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We characterize the expected statistical errors with which the parameters of black-hole binaries can be 
measured from gravitational-wave (GW) observations of their inspiral, merger and ringdown by a network of 
second-generation ground-based GW observatories. We simulate a population of black-hole binaries with uniform 
distribution of component masses in the interval (3, 80) Mq, distributed uniformly in comoving volume, with 
isotropic orientations. From signals producing signal-to-noise ratio > 5 in at least two detectors, we estimate the 
posterior distributions of the binary parameters using the Bayesian parameter estimation code LALInference. 
The GW signals will be redshifted due to the cosmological expansion and we measure only the “redshifted” 
masses. By assuming a cosmology, it is possible to estimate the gravitational masses by inferring the redshift 
from the measured posterior of the luminosity distance. We find that the measurement of the gravitational 
masses will be in general dominated by the error in measuring the luminosity distance. In spite of this, the 
component masses of more than 50% of the population can be measured with accuracy better than ~25% using 
the Advanced LIGO-Virgo network. Additionally, the mass of the final black hole can be measured with median 
accuracy ~ 18%. Spin of the final black hole can be measured with median accuracy ~ 5% (17%) for binaries 
with non-spinning (aligned-spin) black holes. Additional detectors in Japan and India significantly improve the 
accuracy of sky localization, and moderately improve the estimation of luminosity distance, and hence, that of all 
mass parameters. We discuss the implication of these results on the observational evidence of intermediate-mass 
black holes and the estimation of cosmological parameters using GW observations. 


I. INTRODUCTION AND SUMMARY 

While we are nearing the centenary year of the first black 
hole solution in General Relativity (GR) discovered by Karl 
Schwarzschild [1], a number of astronomical observations 
have provided compelling, albeit indirect, evidence of the ex¬ 
istence of astrophysical black holes (see [2] for a review). 
These observations strongly point to the existence of at least 
two populations of black holes: stellar-mass black holes with 
masses ~ 5 - 30 Mq found in x-ray binaries and supermassive 
black holes with masses ~ 10^-10*° Mq found in galactic 
nuclei [2]. There is also suggestive, but much less robust, ev¬ 
idence of a third population of black holes with intermediate 
(~ 10^ - 10^ Mq) masses (see, e.g., [3, 4]). 

The astronomical observations used to infer the existence 
of black holes are indirect. The mass measurements are based 
on observing the motion of luminous objects (such as stars or 
accreting matter) around a compact object. Furthermore, they 
only point to the existence of compact objects which are signif¬ 
icantly more massive than other known compact objects, such 
as neutron stars. These observations are however inadequate to 
establish that these objects are indeed black holes predicted by 
General Relativity as opposed to more exotic compact objects, 
such as Boson stars. Additionally, these (indeed remarkable) 
measurements can be affected by systematic errors that are 
hard to track down, and hence are typically not included in 
the error estimates. Furthermore, the dynamical mass mea¬ 
surements from x-ray binaries only provide lower limits of the 
black hole’s mass since what is measured is the mass function 
of the black hole, which is a combination of the gravitational 
mass and inclination angle [2]. 

The anticipated gravitational-wave (GW) observations by 
the upcoming GW observatories provide a unique opportunity 
to directly measure the masses (and spin angular momenta) 
of black holes in coalescing black-hole binaries. Such binary 
systems are among the prime sources for the first direct de¬ 


tection of GWs using interferometric GW detectors, such as 
Advanced LIGO [5], Advanced Virgo [6], KAGRA [7] and 
LIGO-India [8]. The GWs, produced purely by the motion of 
the black holes in the binary and well described by the GR, 
travel to the detector completely unaffected by the intervening 
matter. By comparing the observational data with theoretical 
templates of the expected signals (as computed by GR) it is pos¬ 
sible to extract the parameters of the binary (such as the masses 
and spin angular momenta of the black holes, sky-location and 
luminosity distance to the binary, etc.). Such observations, 
where the systematic errors are small compared to electro¬ 
magnetic observations, are expected to provide the first direct 
measurements of the mass and spin angular momenta of black 
holes in the next decade. In principle, such observations will 
also provide us an opportunity for testing whether astrophys¬ 
ical black holes behave according to the black hole solutions 
predicted by GR (see, e.g., [9]). 

This paper aims to characterize the limiting statistical er¬ 
rors in the estimation of the parameters of coalescing binary 
black holes by advanced GW detectors. The recent advances 
in numerical relativity [10-12] and analytical relativity [13] 
have provided us with waveform templates that model the com¬ 
plete inspiral, merger and ringdown of the coalescence [14—25]. 
Previous estimates employing Fisher matrix formalism [26] 
have shown that this will significantly improve the parameter 
estimation accuracies of “high-mass” binaries as compared 
to estimates using templates only modeling the inspiral stage 
of the coalescence. Some of the more recent studies employ¬ 
ing Bayesian parameter estimation codes have further demon¬ 
strated this [27-29]. 

Here we present a comprehensive study of the expected sta¬ 
tistical errors in the context of the upcoming advanced GW 
detector networks. We simulate a population of black-hole 
binaries with uniform distribution of component masses in the 
interval (3,80) Mq, distributed uniformly in comoving volume 
with isotropic orientations. The choice of the mass range is 
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motivated by some of the recent population synthesis mod¬ 
els [30]. From signals producing signal-to-noise ratio > 5 in at 
least two detectors, we estimate the posterior distributions of 
the binary parameters using the Bayesian parameter estimation 
code LALInference [31]. In addition to presenting error esti¬ 
mates in the parameters of the binary such as the component 
masses, we also investigate how well the mass and the spin of 
the newly formed black hole (as the product of the merger) can 
be estimated. 

Figures 1 and 4, and Table I summarize the statistical errors 
with which various parameters can be measured from the simu¬ 
lated population of binary black holes. A brief summary of the 
main results is as follows: We find that the error in the mea¬ 
surement of the mass parameters are, in general, dominated 
by the error in the measurement of the luminosity distance to 
the source. This is due to the fact that the masses appear in 
the GW signal completely degenerate with the cosmological 
redshift z, and we can measure only the “redshifted” total mass 
s M(1 H- z). In order to estimate the real gravitational 
mass M, the redshift has to be estimated independently. In 
the absence of any independent measurement, z can be com¬ 
puted from the measured luminosity distance if we assume a 
cosmology. Since luminosity distance is not a very well mea¬ 
sured quantity, the errors in the measurements of the physical 
masses are dominated by the error in measuring the luminosity 
distance. In spite of this, the component masses of more than 
50% of the population of non-spinning can be measured with 
accuracy better than ~25% using the Advanced LIGO-Virgo 
three-detector network. Additionally, the mass and spin of the 
final black hole can be measured with median accuracy ~ 18% 
and ~ 5%. This will be of great interest to astrophysics since 
we are measuring the mass and spin of individual black holes 
(as opposed to combinations like chirp mass or total mass). If 
the black holes have non-precessing spins, the accuracies are 
slightly worse due to the correlation between mass and spin 
(in particular, the estimate of the final spin, which is reduced 
to ~ 17%). 

As far are the redshifted parameters are concerned, as found 
in previous studies [26, 28, 29], we find that the total mass of 
the binary = m\ H- is a better estimated mass parameter 
for heavier binaries (M‘ > ISOMq), while the chirp mass = 
is the best estimated parameter for lighter binaries 
(M^ < 150Mo), where rf = mim 2 lM^ is the symmetric mass 
ratio of the binary. In addition, the mass of the final black 
hole is measured with remarkable accuracy. In fact, we find 
that for masses > 200Mq, the mass of the final black 
hole My is measured to a better precision than even the total 
mass M' . In contrast, the mass of the individual black holes 
is measured rather poorly, limited by the measurement of the 
mass ratio. 

Additional detectors in Japan and India significantly improve 
the accuracy of sky localization, and moderately improve the 
estimation of luminosity distance, and hence that of all mass 
parameters. The accuracy of estimating the physical mass pa¬ 
rameters can be significantly improved by an independent mea¬ 
surement of z, for example from spectroscopic studies of the 
host galaxy, if the host galaxy can be identified. Alternatively, 
an improved measurement of the luminosity distance from GW 
observations (e.g., by reducing the correlation between the 
luminosity distance and the inclination angle by employing 


Parameter 3 detector 5 detector 


Component masses mi ,2 
Chirp mass M 
Total mass M 
Final mass My 
Mass ratio q 
Symmetric mass ratio q 
Final spin Of/Mf 
Luminosity distance di 
Sky location fl (sq. deg.) 


24.5% (41.9%) 
17.4% (29.9%) 
18.2% (30.5%) 
18.3% (31.0%) 
34.8% (60.7%) 
5.4% (16.0%) 
4.5% (13.7%) 
46.6% (75.8%) 
31.25 (121.34) 


22.6% (40.0%) 
15.4% (26.6%) 
15.9% (27.5%) 
16.1% (27.8%) 
33.1% (60.8%) 
5.5% (16.1%) 
4.6% (13.5%) 
40.9% (65.0%) 
8.90 (36.41) 


TABLE I: Summary of the parameter estimation capabilities of the 
upcoming GW detector networks. The table shows the 50 percentile 
(90 percentile) values of Icr errors in estimating the parameters of the 
simulated binary black-hole population using the 3 detector network 
and the 5 detector network considered in the paper. Errors in all the 
parameters except the sky location are fractional errors. Here the 
initial black holes are assumed to have negligible spins. 


waveform templates including the effect of non-quadrupole 
modes, spin precession, etc.) can improve the accuracy of the 
mass parameters. 

The rest of the paper is organized as follows: Sec. II presents 
details of the methodology used in this study, including the 
Bayesian parameter estimation pipeline, gravitational wave¬ 
form models, the astrophysical setup, detector models, etc. 
Section III presents a discussion of the results while Sec. IV dis¬ 
cusses the implication of our results along with some conclud¬ 
ing remarks. Throughout the paper we assume a flat ACDM 
cosmology with parameters Hq - 70km/s/Mpc, fi,„ = 0.3, 
Qa = 0.7, and geometric units: G - c - 1. 


II. METHODOLOGY 

A. Bayesian approach to parameter estimation 

In the presence of a GW signal, the strain of a GW detector 
d{f) is given by 


d(t) ^ n(t) + h(t) (2.1) 


where n(t) is the detector noise time series and h(t) is the GW 
signal. As customary, we are assuming the detector noise to be 
well modelled by a stationary Gaussian process with zero mean. 
Stationarity implies that in the frequency domain all frequency 
components are independent of each other, thus allowing us to 
write the autocorrelation of the noise as 

( 2 . 2 ) 

where we introduced the (one-sided) power spectral density 
S if). Given the above assumptions, the likelihood for a given 
noise realization is given by 


p(n\l) oc exp 


(«|w) 

2 


(2.3) 


where we introduced the scalar product 


ia\b) = 2Re 



a{fW{f) + a*(f)b{f) 

S{f) 


(2.4) 
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FIG. 1: Median values of the Icr errors in estimating various parameters of non-spinning binary black holes using the three-detector (Advanced 
LIGO-Virgo) network as a function of the (redshifted) total mass of the binary. The error in sky location is in square degrees, while the rest are 
fractional errors. The error in estimating the physical mass parameters (left plot) are dominated by the error in estimating the luminosity distance. 
Hence they are estimated with poorer accuracy as compared to their redshifted counterparts. For the case of the redshifted mass parameters 
(middle plot), the chirp mass is the best estimated mass parameter at low masses (M' < 150Mo), while at high masses (M- > 200Mo) the 
total mass M' and the final mass are better estimated. 


The noise statistical properties allow us to write the likelihood 
of a given detector strain realization given a GW signal as 


pid\h, I) oc exp 


id - h\d - h) 
2 


(2.5) 


Assume that the GW signal h{t) depends on some parameters 
9 which are unknown and we intend to infer from the data. We 
can do so using Bayes’ theorem: 


^ ^ p(d\e, h, I) 

pm,h, I) = p(0\h, I) ;/ (2.6) 

p{d\h, I) 

where p0\h, I) is the prior probability distribution for the pa¬ 
rameters 9 given a GW model h, p{d\9, h, I) is the likelihood 
we introduced in Eq. (2.5) and p(d\h, I) is sometimes referred 
to as the evidence or marginal likelihood and it is given by 

p{d\h,l)^ Jd9p(9lh,I)p(dl9,h,I). (2.7) 

Being able to infer the values of the parameters amounts to be¬ 
ing able to construct the posterior distribution p(9ld, h, I). This 
is a formidable task which can only be tackled numerically 
through stochastic samplers. This is mainly due to the large 
dimensionality of the parameter space: For non-spinning bi¬ 
nary systems the vector is 9-dimensional, for more complicated 
(and realistic models) the dimensionality rises quickly to 15 for 
fully spinning stellar systems to even larger integer once more 
detailed physics is included in the gravitational waveform. 

For our purposes, we relied on the LALInference [31] 
stochastic samplers available within the LIGO Algorithm Li¬ 
brary (LAL) [32] . In particular we made use of the LALIn- 
EERENCENEST Software which implements a nested sampling 
algorithm [33] in the context of GW data analysis. Nested 
sampling algorithms aim at solving the integral Eq. (2.7) rather 
than producing samples from Eq. (2.6) and obtain posterior 
samples as a byproduct. We obtain posterior distributions on 


the following parameters this way: the total mass M = nii +m 2 , 
mass ratio q = m\lm 2 , time of arrival at geocenter fo, phase 
of the waveform at a reference frequency location of the 
binary on the sky (a, 6), orientation of the binary w.r.t. the line 
of sight (i, tf/), and the luminosity distance Di. Most of the cal¬ 
culations in this paper are done assuming that the initial black 
holes have negligible spins. However, we do perform one set 
of simulation ascribing non-precessing spins to the initial black 
holes. For this case, the parameter space consists of two addi¬ 
tional dimensions x\ xi corresponding to dimensionless 
spins of the two black holes (aligned/anti-aligned to the orbital 
angular momentum). Orbital eccentricity is assumed to be 
negligible. Since the mass and the spin of the final black hole 
is uniquely predicted by the initial masses and spins, we can 
estimate the posteriors of the mass Mj and dimensionless spin 
afIMf of the final black hole using fitting formulas calibrated 
to numerical-relativity simulations. For this work, we used the 
fitting formulas given in [34] . 


B. Waveform models 


In this paper we characterize the expected statistical errors 
in the estimated parameters by employing the waveform family 
IMRPhenomB [18]. These waveforms describe the GW signals 
from the inspiral, merger and ringdown of binary black holes 
with non-precessing spins. The waveform is written in the 
frequency-domain as h(f) = A(f) where the amplitude 

A(f) and the phase 'V(f) are defined as 


A(f)=Cm 



/o < / < /l 

[l+zh eiv‘] 

/l < / < /2 

U}r£if,f2,0-) 

/2 < / < /3, 


and 


'?(/) = 27r/fo + ^o + 


\2%qiP 




k=2 


( 2 . 8 ) 


(2.9) 
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Above, 


C 


2D, ^ 




- 7/6 


t2/3 


+ cos^ t) D+(a, 5, + 4 cos^, Dx(a, 6 , 



( 2 . 10 ) 


is the amplitude-scaling factor that depends on the antenna pat¬ 
tern functions F+ y^{a, 6 , 41 ) and the inclination angle t, where 
the angles a, 6 describe the sky-location of the binary and [}/ 
the polarization angle. Additionally, /o is the low-frequency 
cutoff of the detector noise where the PSD raises sharply due 
to seismic noise, /a is the high-frequency cutoff above which 
the power in the signal is negligible, /2 and /3 are the transition 
frequencies between the inspiral and the merger and between 
the merger and the ringdown. Also /' s ///j is a dimension¬ 
less frequency, v = 02 = -323/224 + 45 I 77 /I 68 and 

- (27/8 - llril 6 ) are the post-Newtonian (PN) corrections 
(1.5PN accurate) to the leading order amplitude of the inspiral, 
ei = 1.4547;^^ - 1.8897 and 62 = -1.8153;\/ H- 1.6557 are pa¬ 
rameters describing the amplitude of the merger obtained from 
numerical-relativity simulations. ■C{f,f2,cr) is a Lorentzian 
function with width cr centered around /2 while Wm and Wy 
makes the amplimde continuous over the transition frequen¬ 
cies /i and / 2 . In the definition of the phase, fo is the arrival 
time of the signal at the detector, i/jq is the phase at a reference 
frequency and xpk are phenomenological parameters describing 
the phase evolution of the binary. The phenomenological pa¬ 
rameters {/i, / 2 , /s, cr, i/,j.) are entirely functions of the physical 
parameters [M, t],x] and are given in Table I of [18]. Note that 
in this section, we did not include the effect of cosmological 
redshift in the observed waveform. This is described in the 
next section. 


C. Astrophysical set up 


Our purpose is to understand the accuracy with which we can 
realistically expect to measure the parameters of binary black 
holes from GW observations of high-mass stellar black holes. 
For this purpose we simulate a population of 5000 binaries 
with rest-frame component masses uniformly distributed in 
(3,80) Mq. The binaries are also distributed uniformly in co¬ 
moving volume, therefore uniform on the celestial 2-sphere 
and in redshift according to the probability distribution [35] 


dR(z) 1 
dz RiZm-ix) 


( 2 . 11 ) 


in the redshift z space. In Eq. (2.11), R{z) is the cosmic rate 
whose evolution with redshift is equal to 


dR(z) _ dV roe(z) 
dz dz 1 H- z 


( 2 . 12 ) 


where rg is the rate of binary-black-hole coalescence in the 
local Universe, e(z) is the cosmic star formation rate and V 
is the co-moving volume. In a Friedmann-Robertson-Walker- 
Lemaitre universe 



Dliz) 

(l+z)^H(z) 


(2.13) 



FIG. 2: The top plot shows the horizon distance (distance to which 
optimally located- and oriented binaries can be observed with an 
optimal signal-to-noise ratio of 8) of the Advanced LIGO detector 
towards non-spinning, equal-mass binaries binaries as a function of 
the (redshifted) total mass of the binary. The bottom plot shows the 
expected distribution of SNRs from the population of binary black 
holes that cross the detection criterion. 



where 


H{z) = Ho VDm(l +z)3 -hQ^(1 -hz)2-hQa (2.14) 

is the Flubble parameter. Above, Hq is the Hubble constant, 
Q,„ is the fractional matter density, is the fractional energy 
density due to a cosmological constant, = 1 - Q„, - Qa is 
the fractional curvature energy density, which we take to be 
zero (Q^ = 0). The luminosity distance is then given by [36] 

dz' 

Dl(z) = (1+z) . (2.15) 

Jo H{z') 

Since we are not interested in predictions relative to the coales¬ 
cence rates of the high mass systems, we considered here the 
case ro - 1 and e(z) = const. 

It is worth noting that the mass parameters and the red¬ 
shift are in general completely degenerate [37], thus what is 
observed in the detector is not the rest-frame mass, but the “red- 
shifted” one. Therefore, the mass distribution of all detected 
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mi,2 [Mq] m\ 2 [Mq] M' [Mq] q 


FIG. 3: The left plot shows the distribution of the rest-frame component masses m, 2 of the simulated population of non-spinning binaries 
(light colored histograms) and that of the sub-population that crossed the detection threshold for the three-detector network (dark colored 
histograms). The two middle plots show the distributions of the redshifted component masses m\ ^ and the total mass M-. Note that, even 
stellar-mass (mj 2 < 80Mq) black holes can appear to have intermediate masses {m\ ^ > lOOMg) due to the degeneracy between mass and 
cosmological redshift. The right plot shows the distributions of the mass ratio q = m 2 1 mi. Note that comparable-mass {q » 0) binaries are 
detected preferentially due their larger intrinsic luminosity. 


sources will be quite different from the intrinsic distribution 
(see Fig. 3). We will denote the redshifted mass parameters 
with a superscript z {m\ 2 , etc.). 

For the same set of masses, sky position, distances, orienta¬ 
tion angles, we considered two scenarios, one where the spins 
of the BFIs are negligible and one where the dimensionless 
spins can assume take arbitrary values in the interval (-1,1), 
but are restricted to lie in a direction aligned/anti-aligned to 
the orbital angular momentum, so that the spins do not pre- 
cess. We did not consider the general case of fully precessing 
binary systems as no fully precessing waveform including also 
merger and ringdown is available yet and, considering the 
mass range we are interested in, merger and ringdown con¬ 
tribute significantly if not predominantly to the accumulation 
of signal-to-noise ratio in the detector. 


D. Detector models 


We consider two upcoming/anticipated networks of ad¬ 
vanced GW observatories: 1) a three-detector network con¬ 
sisting of two Advanced LIGO detectors located in Hanford, 
WA and Livingston, LA in the United States and one Advanced 
Virgo network located in Cascina, Italy. 2) a five-detector net¬ 
work consisting of the three Advanced LIGO-Virgo detectors 
mentioned above and the KAGRA detector located in Kamioka, 
Japan and the LIGO-India detector located in India. Since the 
location of the LIGO-India is not finalized, we assume a fidu¬ 
cial location north of Bangalore [38]. The Advanced Virgo 
detector is assumed to have the sensitivity given in [39] while 
the rest of the detectors are assumed to be having the sensitiv¬ 
ity of the Advanced LIGO in the “high-power-zero-detuning” 
configuration [40]. We assume a low-frequency cutoff of 10 
Hz below which the detector has extremely poor sensitivity 
due to seismic noise. 


E. Simulations 

From the simulated population of binary black holes de¬ 
scribed in Sec. IIC, we select those binaries producing signal- 
to-noise ratio > 5 in at least two detectors. This is to simulate 
a realistic detection scenario with two-detector coincidence. 
The distribution of the network signal-to-noise ratio of these 
“detected” binaries is shown in the bottom panel of Fig. 2. The 
top panel of this figure shows the horizon distance (the distance 
at which the detector can observe an optimally-located- and 
oriented binary with signal-to-noise ratio of 8) of Advanced 
LIGO towards non-spinning equal-mass binaries, which gives 
an idea of the distance reach of the detector network. Posterior 
distributions of all the binary parameters are computed using 
the Bayesian parameter estimation code LALInferencenest. 
We choose prior distributions as follows: observed component 
masses uniform in the interval (2,300) Mq, sky location and 
orientation uniform on the 2-sphere, a uniform distribution for 
the time of arrival at geocenter to with a width of 0.1 s around 
the true arrival time, phase cpo uniformly distributed between 0 
and Itt and for the luminosity distance we assume a prior distri¬ 
bution given by Eqs. (2.11), (2.12), (2.13) withD/, 6 (1,50000) 
Mpc. For the spinning case, we assume a uniform prior on 
the dimensionless spins for both components between [-1,1]. 
From the marginalized posterior of each parameter of interest, 
we compute the Icr conhdence intervals. 


III. RESULTS AND DISCUSSION 

The main results of our work are summarized in Figs. 1 
and 4, and in Table 1. To remind the reader, we simulate an 
astrophysical simulation of 5000 black holes with component 
masses uniformly distributed between (3,80) Mq. Figure 1 
shows the median value of the Icr errors on various estimated 
parameters as a function of the (redshifted) total mass M' of 
the binary. Consistent with our expectations, we find that the 
chirp mass is the most precisely estimated quantity for low 
masses < 150Mo, the total mass M' best estimated for 
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AM'/M' AM^IM- AM^My 


FIG. 4: The top panel shows the distribution of the Icr errors in estimating the “redshifted” mass parameters (redshifted component masses m\ 
chirp mass M', total mass M' and the final mass My) from the simulated population of coalescing binary black holes. In each subplot, the three 
distributions correspond to non-spinning binary black holes detected by the 3-detector network and the 5-detector network that we consider, 
and the population of binary black holes with non-precessing spins detected by the 3-detector network. The vertical lines show the median of 
the distribution. The bottom panels show the distribution of the errors in estimating the “physical” parameters. It can be seen that the error in 
measuring the “physical” mass parameters AI, M and Mf are significantly larger than that in measuring the “redshifted” parameters, and are 
dominated by the error in measuring the luminosity distance. The error in the component masses are dominated by that in measuring the mass 
ratio of the system, and hence are not very different in the top and bottom panels. Note that the mass ratio rj, the final spin Uf/Mf, the luminosity 
distance di^ and the sky location are independent of the cosmological redshift. 






Ami_2/n2i,2 


AM/M 
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sky loc. (sq. deg.) Adi/di Arj/r] Aaf/of 


intermediate masses 150/Wo < < 200Mq, while the mass 

of the final black hole /W^. is the best estimated mass parameter 
for high masses > 200/Wo. Figure 4 shows the distributions 
of the expected Icr statistical errors on the estimated quantities 
from the simulated source population. The top panel shows 
the distributions of the “redshifted” parameters and the bottom 
panel shows the distribution of the “physical” parameters. The 
vertical lines of the same colors indicate the median of the 
respective distributions. A verbal summary of the results is as 
follows; 

a. Measurement of redshifted parameters: For the case 
of non-spinning binaries, the chirp mass the total mass 
/W^, and the final mass /W^ are estimated to better than 4% 


(4%) for more than half of the population, using a five (three) 
detector network. In contrast, the component masses oti _2 
are estimated to about 16% (17%). In the case of aligned- 
spin binaries, the measurement accuracies deteriorate, due to 
correlations between mass and spin. The median measurement 
accuracies for m\ Ai^, M^, My are 19%, 9%, 9%, 8% using 
a three-detector network. 

b. Measurement of “physical” parameters: For the case 
of non-spinning binaries, the chirp mass A4, the total mass 
M, and the final mass Mf are estimated to better than 16% 
(18%) for more than half of the population, using a five (three) 
detector network. Component masses mi 2 are estimated to 
about 23% (25%). For aligned-spin binaries, the median mea- 
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FIG. 5: Typical posteriors obtained using LALInference. The solid red lines indicate the injection values. The dashed magenta, cyan and yellow 
contours indicate the 68.3%, 95.5% and 99.7% confidence regions respectively. The left panel shows the correlation between the component 
masses. For large masses, the total mass is roughly the most precisely estimated quantity. The posterior samples are thus correlated along 
the M~ = constant line. The right panel shows the correlation between luminosity distance and inclination. In this case di/jl + cos^ t) is the 
precisely estimated quantity. Because of this degeneracy, neither of the two physical quantities are well-estimated. 


surement accuracies for mj 2 , M., M, Mf are 25%, 18%, 19%, 
20% using a three-detector network. 

The sky location, as expected, shows a significant improve¬ 
ment with 5 detectors as compared to 3 detectors. This is 
because the leading estimate of the sky location comes from 
the difference in the time of arrival of the signal at the different 
detectors, and with 3 detectors, in many cases, one cannot dis¬ 
tinguish between two antipodal points in the sky. The estimate 
on the luminosity distance is quite imprecise as compared to 
the redshifted masses above — the median error on distance 
is about 40% and 50% with 5 and 3 detectors respectively. 
Since the distance is used to calculate the redshift in order to 
translate from the observed to the physical masses, the errors 
in the physical chirp, total, and final masses are dominated by 
the inaccuracy of the distance estimate, and are much worse 
compared to their corresponding redshifted values. On the 
other hand, the mass ratio is unaffected by redshift. The errors 
in the component masses do not change much between their 
redshifted and physical values, since they are dominated by the 
error in estimating the mass ratio. Moreover, since the final 
spin afIMf is determined entirely by the mass ratio and the 
initial (dimensionless) spins, the final spin is also unaffected 
by redshift. 

We see that allowing for the possibility of aligned-spins, in 
general, deteriorates the estimation of mass and spin parame¬ 
ters. This is due to the correlation between the mass and the 
spins. This degeneracy is particularly pernicious for the case of 
templates with aligned spins [41], which could be improved by 
including the effect of spin precession [42]. The measurement 
of the physical mass masses, M,M, and Mf can be improved 
by improving the measurement of luminosity distance. Lumi¬ 
nosity distance is not a well-estimated quantity, mainly because 
of its correlation with the inclination angle (see Figure 5). The 
degeneracy can be slightly broken and the estimates made more 


precise by including more detectors. We already see that going 
from 3 to 5 detectors improves the distance estimate by a factor 
of about 1.5. Including higher harmonics in the waveform can 
be even more effective in breaking this degeneracy, and can be 
of great importance in the measurement of the physical masses 
of black holes. 

Let us also list some of the caveats of our results; Note that 
we used the analytical IMR gravitational-waveform family IM- 
RPhenomB that has been calibrated against numerical-relativity 
simulations with mass ratio q > 1/4. However, the effect of 
any waveform systematics is likely to be very small: As seen 
in the right plot of Fig. 3, the population of detected binaries 
is significantly dominated by those with comparable masses 
(due their larger intrinsic luminosity). Additionally, in our 
simulation we use a uniform distribution of component spins 
in the interval (-1,1), while IMRPhenomB has been calibrated 
against numerical-relativity waveforms with moderate spins 
(\Xi, 2 \ ^ 0.85 for q' = 1 and \xi, 2 \ ^ 0.75 for q < 1). This 
can introduce non-negligible systematics in the case of highly 
spinning binaries. Even though it is unlikely that this will 
significantly change the median values of the statistical errors 
that we quote, the error estimates for spinning binaries need to 
be looked at keeping this caveat in mind. 


IV. IMPLICATIONS OF THE RESULTS 

Measuring black-hole mass- and spin functions :— We find 
that the errors in measuring the “physical” mass parameters 
is in generally dominated by the error in measuring the lumi¬ 
nosity distance, which is often overlooked. Even when taking 
this effect into account, we see that it is possible to measure 
the component masses of more than ~ 50% of the assumed 
population with accuracies better than 25%. We also find that 
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the mass of the final black hole can be measured with a slightly 
better median accuracy of ~ 18%. More interestingly, the 
spin of the final black hole can be measured with a remark¬ 
able accuracy of ~ 5 - 17%. This has to be combined with 
the fact that GW-based measurements will have significantly 
smaller systematic errors as compared to the electromagnetic 
measurements. 

Detecting intermediate-mass black holes: — One of the in¬ 
teresting scientific potential of GW observations of binary 
black holes is to shed light on the possible existence of 
intermediate-mass black holes [28, 29]. Note that, due to 
cosmological redshift, even stellar-mass black holes (mi ,2 ^ 
8 OM 0 ) can have observed mass greater than IOOMq (see Fig. 3). 
Hence the mass measurements from GW observations have to 
be interpreted carefully. However, it is possible to measure the 
physical component masses with reasonable accuracy (~ 25%). 
More interestingly, since we are able to measure the mass of 
the final black hole well, GW observations might enable us 
to witness the birth of intermediate-mass black holes through 
mergers. This would help shed light on the the yet unclear 
channels of formation. 

Measuring cosmological parameters: — It is, in principle, 
possible to estimate some of the cosmological parameters (such 
as Hq and Q^) by combining the estimates of the luminosity 
distance with independent estimates of the cosmological red- 
shift from spectroscopic observations of the host galaxy [43]. 
Using the median error estimates of the luminosity distance 
and sky location (~ 50% and 30 sq. deg for the three-detector 


network), we estimate that ~ 200 potential host galaxies are 
likely to be present in the allowed confidence region. This 
will be considerably reduced in the case of a five-detector net¬ 
work. Still, it will be impossible to uniquely identify the host 
galaxy. However, reasonable constraints on the cosmological 
parameters can be obtained by combining multiple GW obser¬ 
vations [44]. This is being explored in detail in an ongoing 
work [45]. 
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